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We investigate the power of quantum systems for the simulation of Hamiltonian time evolutions 
on a cubic lattice under the constraint of translational invariance. Given a set of translationally 
invariant local Hamiltonians and short range interactions we determine time evolutions which can 
and those that can not be simulated. Whereas for general spin systems no finite universal set of 
generating interactions is shown to exist, universality turns out to be generic for quadratic bosonic 
and fermionic nearest-neighbor interactions when supplemented by all translationally invariant on- 
site Hamiltonians. 



I. INTRODUCTION 

One of the most promising applications of quantum 
computers is that of being a working horse for physicists 
who want to determine the time evolution of a theoreti- 
cally modeled quantum system. On a classical computer 
this is a daunting task burdened by the notorious ex- 
ponential growth of the underlying Hilbert space. As 
pointed out by Feynman [l], quantum systems can, 
however, be used to efficiently simulate quantum mechan- 
ical time evolutions — provided that we have sufficient co- 
herent control on the system. In this direction enormous 
progress has been made during the last years, in partic- 
ular in systems of optical lattices [H and ion traps M. 
Moreover, it was realized that quantum simulators [6| are 
much less demanding than quantum computers and, in 
fact, pioneering experiments simulating quantum phase 
transitions in systems of cold atomic gases [3| have al- 
ready turned some of the visions Q into reality. 

One of the fundamental questions in the field of quan- 
tum simulations is the following: Given a set of inter- 
actions we can engineer with a particular system, which 
are the Hamiltonians that can be simulated? Concerning 
gates, i.e., discrete time unitary evolutions, it has been 
shown in the early days of quantum information theory 
that almost any two-qubit gate is universal [8| . Similarly, 
any fixed entangling two-body interaction was shown to 
be capable of simulating any other two-body Hamiltonian 
when supplemented by the set of all local unitaries [9(. 
The many-body analogue of this problem was solved in 
(To| and the efficiency of quantum simulations was stud- 
ied in various contexts (cf. fill]). 

All these schemes are based on the addressing of sites, 
i.e., local control. Imagine now that we have a chain in 
which we cannot address each particle individually but 
only apply global single-particle and nearest-neighbor in- 
teractions. Can we simulate the evolution of a next-to- 
nearest neighbor interaction Hamiltonian, or obtain some 
long-range (e.g., dipole) coupling, or even a three-particle 
interaction Hamiltonian? 

In this article we will concentrate on the case in which 
the interactions at hand are short range and translation- 
ally invariant as it is (approximately) the case in different 
experimental set-ups, like in the case of atoms in optical 
lattices or in many other systems that naturally appear in 



the context of condensed matter and statistical physics. 
In order to make the problem mathematically tractable 
and to exploit its symmetries we will consider periodic 
boundary conditions, even though typically physical sys- 
tems have open ones. In this sense, our results may not 
be directly applicable to certain physical situations. In 
any case, we expect that our work will be a step forward 
to the establishment of what can and cannot be simu- 
lated with certain quantum systems. We will consider 
three different systems: spins, fermions and bosons. A 
summary of the results of this work is given in the fol- 
lowing section. 

II. STATEMENT OF THE PROBLEM AND 
SUMMARY OF RESULTS 

Consider a cubic lattice of N sites with periodic bound- 
ary conditions in arbitrary spatial dimension. Assume 
that we can implement every Hamiltonian from a given 
set iS = {7^1,7^2, ■ ■ ■} of translationally invariant Hamil- 
tonians and in this way achieve every unitary time evolu- 
tion of the form e lHjt for arbitrary t € R. Note that this 
assumes that both ±Tlj are available. The question we 
are going to address is, which evolutions can be simulated 
by concatenating evolutions generated by the elements of 
S. Our main interest lies in sets which contain all on-site 
Hamiltonians and specific nearest-neighbor interactions. 

The natural language for tackling this problem is the 
one of Lie algebras [H, [l3| since the set of reachable 
interactions is given by the Lie algebra C generated by 
the set iS. This follows from the Lie -Trotter formulae 

m 

e aL k +PL t = Hm / a L t /n e ^,/„r a, [3 £ R, 
n — >-oo V / 

e \LuM\ = lim f e L k /V^ e Li/V^ e -Lk/V^ e -L,/V^Y i 
n — >oo V / 

where is a representation of the generator iHk- When 
applying the Lie- Trotter formulae to the elements of iS 
we can obtain all commutators and real linear combina- 
tions of its elements, i.e., we end up with the Lie algebra 
generated by S. Conversely, it follows from the Baker- 
Campbell-Haussdorff formula [l5[ that all simulatable in- 
teractions can be written in this way. We will study the 
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cases of £>-dimensional 'spin' systems (C C su d n) as 
well as quadratic Hamiltonians in fermionic (£ C sc>2n) 
and bosonic (£ C sp2n) operators. The following gives 
a simplified summary of the main results. Hamiltonians 
and interactions are meant to be translationally invari- 
ant throughout and it is assumed that interactions along 
different directions can be implemented independently. 

• Fermions: All simulated evolutions have real tun- 
neling/hopping amplitudes. Within this set generic 
nearest-neighbor interactions are universal for the 
simulation of any translationally invariant interac- 
tion when supplemented with all on-site Hamilto- 
nians. Whereas for cubes with odd edge length the 
proof of universality requires interactions along all 
axes and diagonals, the diagonals are not required 
for even edge length. 

• Bosons: All simulated evolutions are point sym- 
metric. Within this set every nearest-neighbor in- 
teraction (available along axes and diagonals) is 
universal for the simulation of any translationally 
invariant interaction when supplemented with all 
on-site Hamiltonians. 

• Spins: There is no universal set of nearest-neighbor 
interactions. Moreover, if / is a factor of the edge 
length m of the cubic lattice then there is no univer- 
sal set of interactions with interaction range smaller 
than /. In particular, if m is even, not all next- 
to-nearest neighbor interactions can be simulated 
from nearest neighbor ones. Sets of Hamiltonians 
that can be simulated are constructed. 

Whereas in the case of quadratic bosonic and fermionic 
Hamiltonians a rather exhaustive characterization of sim- 
ulatable time evolutions is possible, a full characteriza- 
tion of simulatable spins systems still remains an open 
problem. 

We will start with introducing some preliminaries on 
quadratic Hamiltonians in Sec lIIIl Scc lIVI will then treat 
fermionic and Sec|V] bosonic systems. Both start with 
the one-dimensional case which is then generalized to ar- 
bitrary dimensional cubic lattices. Finally spin systems 
are addressed in SeclVIl 

III. QUADRATIC HAMILTONIANS 

This section will introduce the basic notions and the 
notation used in Sees. IIVIVI The presentation is a col- 
lection of tools widely used in the literature on trans- 
lationally invariant quasi-free fermionic [HI, [T3, EH fl9j 
and bosonic [lliHH systems. We consider a system of N 
fermionic or bosonic modes characterized by a quadratic 
Hamiltonian 

N 

H = ^ A H a k ai + B k ia k a\ + C k ia\ai + D k ia\a\. (1) 

k,l = l 



Here, a k and a k are creation and annihilation operators 
satisfying the canonical (anti)-commutation relations 

CAR: {a k ,ai} = 0, {a k ,aj} = 5 k i (fermions), (2) 
CCR: [a k ,ai] = 0, [ak,aj]=6ki (bosons). (3) 

By defining a vector a = (ai, . . . , on, a\, . . . aL) and a 
Hamiltonian matrix 

Eq.([T]) can be written in the compact form Ti = aHa T . 
The Hermiticity of TL implies the relations 

B = B\ C = & ,A = Dt. (5) 

We will identify Hamiltonians which differ by multiples 
of the identity as they give rise to undistinguishable time 
evolutions. The commutation relations can then be ex- 
ploited to symmetrize the Hamiltonian matrix H such 
that 

A = tA t , D = tD t , B = tC t , (6) 

where r = 1 for bosons and r = — 1 in the case of 
fermions. Instead of working with 2N creation and anni- 
hilation operators it is often convenient to introduce 2N 
hermitian operators c k via 

Cfc = (<4 + a k)/ V%, c k+N = i{a\ - a k )/ y/2. (7) 

In the case of fermions these are the Majorana operators 
obeying the anti-commutation relation 

{c k ,c t }=6 kl . (8) 

For bosons the c k are the position and c k+ N momen- 
tum operators, and the commutation relations can be 
expressed in terms of the symplectic matrix a via 

[c k , Cl ] = ia kl , o- = ( ° t I ) , 1 e R NxN . (9) 

Eq. ([1]) can now be written in the form 

k,l ^ ' 

Exploiting again the commutation relations we can 
choose the Hamiltonian matrix H real and (anti-) sym- 
metric with H = tH t . The Hamiltonian matrices of the 
two representations are related via 

fr _ Vt ( X-Y-i(W + tW t ) X + Y + i(W - tW t ) \ 
4 \ X + Y - i(W - tW t ) X - Y + i(W + rW T ) )' 

Time- evolution: We are interested in time-evolutions 
generated by quadratic Hamiltonians of the form in Eq. 
(TT|). These are canonical transformations which preserve 



3 



the (anti-) commutation relations and act (in the Heisen- 
berg picture) linearly on the c^'s: 



e mt c k e 



mt 



N 



= > T lk ci 



1=1 



(11) 



In the fermionic case the CAR are preserved iff T E 
0(2N) is an element of the orthogonal group in 2N di- 
mensions. This group has two components corresponding 
to elements with determinant ±1. As time evolution has 
to be in the part connected to the identity (for t = 0) 
we have that T € SO(2N) is an element of the spe- 
cial orthogonal group. For bosons the preservation of 
the commutation relations implies that T is a symplec- 
tic matrix, i.e. ToT T = a. Both groups SO(2N) and 
Sp(2N) are Lie groups and we can express T in terms of 
the exponential map acting on the respective Lie algebra, 
i.e., T — e tL . From the infinitesimal version of Eq. lfTTj) 
we obtain a simple relation between the generator L and 
the Hamiltonian matrix H: 



L = —H for fcrmions 
L = Ha T for bosons. 



(12) 
(13) 



Translational invariant systems: We will throughout 
consider translationally invariant systems on cubic lat- 
tices in d spatial dimensions with periodic boundary con- 
ditions. Hence, the indices of the Hamiltonian matrix 
Hki which correspond to two points on the lattice are d- 
component vectors k, I € where m is the edge length 
of the cube, i.e., N = m d . The translational invariance is 
expressed by the fact that the matrix elements G k i , of the 
blocks G G {X, Y, W} of H depend only on the relative 
distance k — I. Taking into account the periodic bound- 
ary conditions, k — I is understood modulo m in each 
component. Such matrices are called circulant, and we 
will denote by Ca and Cs the set of circulant symmetric 
and antisymmetric matrices, respectively. All circulant 
matrices can be diagonalized simultaneously by Fourier 
transformation 



G = T® d GP® d = diag 



k&<L 



(14) 



1 



T vq = -^e— pq , p,q£Zm, (15) 



where G k = G k ,o is the entry of the A:-th off-diagonal 
of the matrix G. It follows from (TTJJ that all circulant 
matrices mutually commute. 



IV. SIMULATIONS IN FERMIONIC SYSTEMS 

In this section we study the set of interactions that can 
be simulated in a translationally invariant fermionic sys- 
tem starting with quadratic local transformations and 
nearest neighbor-interactions. Making use of the fact 



that the blocks X, Y and W in Eq. (Tl0[) mutually com- 
mute we calculate the commutator L" = [L, L'} of two 
generators L and L' given by 

X( ' ] W( ' ] \ na\ 



-.0 - 



and obtain 



L = 



X" 
-W" T 



W" 

-X" 



X" 
W" 



W'W T - WW' T 
W{Y' - X') (17) 
—W'(Y — X) 

Note that by Eq. pT]) every commutator has the sym- 
metry Y" = —X". Hence, if we start with a set S 
of Hamiltonians with corresponding generators Sc = 
{Li, L2, . . .}, then every element of the generated Lie al- 
gebra C has this form up to linear combinations of el- 
ements in Sc- On the level of Hamiltonians this sym- 
metry corresponds to real tunneling/hopping coefficients 
Bu = -C lk £ M in Eq.Q. We will denote by K the 
vector space of all matrices of the form (|16[) for which 
Y = -X. 

Let us now introduce the elements of the set Sc corre- 
sponding to all local Hamiltonians and specific nearest- 
neighbor interactions. Every generator L of a local 
Hamiltonian is proportional to 



E 



1 

-1 



(18) 



For giving an explicit form to the nearest-neighbor inter- 
action, we define a matrix via 



M. 



(") _ x 

kl ~ °l,k+v, 



(19) 



where v,k,l € and the addition is modulo m in each 
of the components. This leads to the properties 

M^MW = M {V1+V2 \ M<°) = 1. (20) 

Moreover, we define the matrices 

= M (v) + M { ~ v \ = - M ( -- v \(21) 



H 



(v) 
X 



M ( _ v) 
—M 



(v) , 



H 



(v) 



M, 



(v) 



w(±) 



(±) 



M (v)T 
IV1 (±) " 



(22) 
(23) 



where the indices X and W refer to a non-vanishing X- 
and VF-block respectively. Denoting by € TL^ the ba- 
sis vectors = Sij, every Hamiltonian matrix corre- 
sponding to a nearest-neighbor interaction along e, is of 
the form 



if - I X ° W o 
Ho = ' -W? Y 



(24) 



xM. 



(e<) 



(e 4 ) 



where x,y,w,w G R. We will now start studying 
one-dimensional systems and then generalize to the d- 
dimensional case. 
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A. Simulations in one-dimensional fermionic 
systems 

In this section we consider quadratic fermionic Hamil- 
tonians with translational symmetry on a ring of m sites. 
We will give an exhaustive characterization of nearest- 
neighbor Hamiltonians which are universal for the simu- 
lation of all interactions obeying the symmetry Y = —X, 
when supplemented by all on-site Hamiltonians. The re- 
sults depend on whether m is even or odd. 

Theorem 1 Consider a translationally invariant 
fermionic systems of m sites on a ring with periodic 
boundary conditions. Starting with all one-particle 
transformations which are proportional to the matrix 
E defined in (j!8[) and one nearest-neighbor interaction 
Hq of the form (I24[) we can simulate the following set 
of interactions depending on the symmetry properties of 
Xq, Yo and Wq: 

1. m — 2n + 1 odd: 



If Wq $l Cs or Xq ^ Yq, we will show in the first step 
by induction over k that the set X defined in (|2b[) can 

be simulated. For k = 1, we can get H% and H&_ 
by taking the commutator of Hq with the one-particle 
transformation E: If W € C s , then [H ,E] ~ H^ 1 } 
and ffjf 1 can be obtained using (j2"6"l) . If Xq — Yq, then 

[H ,E] ~ H^ l] an 
and Xq ^ Yq, then 



[H Q ,E] ~ H { x l] and we get H^ 1 } by If W $ C s 



[H , E]/(w -w) + [[H , E],E]/2{x Q - y ) ~ 

and according to (|26[) we also get H^K From (|28[) we 
see that we get 

H%f = [HW,Hfrl]/2 + 2E. 
Now let k > 1. Using (J2SJ) and ([22]) we get 

ru-Cei) rr(2fee 1 ) 1 „((2fc-l)ei) _ „((2(fe+l)-l) ei ) 



(a) Xq = Y"o, Wo S Cs.' No further interaction can 
be simulated. 

(b) X ^ F or W £ C s : TTie space ft (i.e. X = 
— Y) can be simulated. 

2. m = 2n even: 

(a) Xq = Y"o; W) £ C5: No further interaction can 
be simulated. 

(b) Wo € or Wq = 0: The space spanned by 



J= {if 



((2fc-l)ei) H (2k ei ) fJ .((2fe-l)ei) 



A' 



-Hi 



'} 



can fee simulated. (For the definition of Hx 
and H w see Eqs. (|2"2")l . (12"5)) ). 

(cj Wo ^ Ca, Hq G TZ: The space ft can fee simu- 
lated. 

(d) Wq 4- Ca, Hq £ ft: TTie space spanned by 



J =\ tl X ■> ti W '% _ / fc 6 N 



which implies that we also get if^ 2 ( fe+1 ' 1 ) ei )_ As 

rrr(ei) rr((2(fc+l)-l)ei)i rr(2ei) _ rr(2(fe+l)ei) 
\n x ,n w _ \^n w+ — n w+ , 

we have shown that we can simulate T. Using the re- 
lations (|2"5l) - (j2"9"|) , we see that I is closed under the 
commutator bracket. 

If m = 2n + 1, I is a basis of all possible interactions 
of the space ft because of the periodic boundary condi- 
tions. To see this, define for an arbitrary k the num- 
ber k' = k + n + 1. As 2fc'modm = 2k + 1, we see 



that H 



((2fc'-l)ei) 



„(2kei) „((2fe'-l) ei ) 
A J-^W- 



rr(2fe' ei ) _ „(2(fe+l)ei) 



H^ l} and 



[ W+ — , which proves (pTbj) . 

Now let m = 2n. If W € or W = 0, then 

• -*+**». a-(-f J M )fi- 

The elements of T are the only ones that can be simulated 



can be simulated. 
Proof For the proof we will need the relations 



[H { x ei) ,E] 


orr( k i) 

— ) 


(25) 


[E,H$ ei) ] 


- ^A > 


(26) 


rrr(fcei) rr(£ei)i 
l H X ) H X I 


= o, 


(27) 


rrr(fcei) rr(£ei)i 
l H X ) J 


- A H w - H w )> 


(28) 


rrr(fcei) rj-(lei)i 




(29) 



For X = Y ,W eC s , [H a , E] = according to (TTTJ) so 
that we cannot simulate any further interaction (up to 
multiples of E). This proves (jTaJ) and (f2"a[) . 



[ff.fffe] - (x - yo )[H^',H^']ex 

where I = 2k, 2k — 1 respectively and H commutes with 



(ei) irfieiji 



((2fe-l)ei) 
A 



This proves (|2b 



If Wo ^ and iJ € ft, then 

H = x H^ - w H&l + (w + wo)H<& ] 
so that we can extract H^\ According to (|28| 



rH .((2fe-l)ei) H -(e 1 ) 1 



0(T j(2ke 1 ) jj 



(_2(fc-l) ei ). 



so that we can get H^ kei ' > as H^) — E, and we can 



simulate H^ kei ^ using (|2"6"|) . It remains to show that we 



5 



can simulate H, 



((2fc-l)ei) 

w 



Note that the possibility of 



simulating H^ ei ^ implies that we can get H^ kei \ as 

H^ + [[H^\e},E]/2 = H^\ 
According to ([25]) 

ro-Cei) Tr(2ke 1 ) ] _ „((2fc+l)ei) „(-(2fc-l) ei ) 
r ((2fc-l)ei) o-(-ei) 



as ii^ is available. This 



so that we can get iJ w 
proves (j2c|) . 

Finally we consider the case where Wq ^ and -ffo ^ 72. 
Then 



■■• • « 



ffo = -ff + ^o-ff^J, 



# = 



We will now calculate the commutator of H with all ele- 
ments of I in order to see if we get additional interactions. 
From 



[H,H X 



« 2k - 1 ^]=-2(w +w )H% k _ ei) 



we see that we can get H^ kei \ as iT^jf G J and using 
(HI we get H ( x kei) . As 



((2/+fe) ei ) o-((2i-fe) ei ) 



we have shown that the set J can be simulated. Using 
(l25f - (l29l) . we see that {J, H} is closed under the com- 
mutator bracket, which proves (|2d[) . □ 



B. Simulations in <i-dimensional fermionic systems 

This section will generalize the previous results to sys- 
tems in d spatial dimensions. The following theorem 
shows that certain nearest-neighbor interactions are uni- 
versal for simulating the space TZ (i.e. Y = —X) on a 
c?-dimcnsional cube. 

Theorem 2 Consider a fermionic systems on a d- 
dimensional translationally invariant cubic lattice with 
m d sites and periodic boundary conditions. Then the fol- 
lowing sets of nearest-neighbor interactions together with 
all on-site transformations are complete for simulating 
the space TZ: 

1. m — 2n + 1 odd: 



FIG. 1: Two-dimensional fermionic lattice with box B2 (see 
proof of Thm. |2j. The interactions along ei and ei (solid 
arrows) are used to simulate the interactions in the direction 
2(ei + €2) (dashed arrow). 



Proof We start with an odd number of fermions, m = 
2n + 1. For the proof we will consider interactions with 
a maximal interaction range in each direction of the 
lattice. To do so, we define for every integer z G N 
the d-dimensional cube of edge length 2z, B z ~ {v G 
Z£J IMU < z}. Then a Hamiltonian fl"^ 00 ) where 



(v) 
X 



ujj ) G B z couples a given lattice site s 
only with sites which lie in a cube of edge size 2z with s 
in its center. We will show by induction over z that H 

and Hyfi can be simulated. We start with a minimal 
edge length of 2, i.e. z = 1 and define Af v — {i | \vi\ = 1} 

with cardinality \M V \- We will show that H x , 
can be simulated for \Af v (i) \ — 1, . . . , d, i.e., for an arbi- 
trary number of non- vanishing components of the vector 
p or |A4(i)| = 1, the vector t>W has only one non- 
vanishing component — ±ej, and the situation is 
as the one of theorem [T] Hence H x and can 
be simulated for arbitrary i. Now let |A/„(i)| = r > 1, 
j G A^ti), i.e. we want to simulate an interaction in 
the direction of the diagonals as depicted in figure Q] As 
\J\f (1) _ (l) I = r — 1 we know by induction over the cardi- 

nality of N v m that H w 3 



we get H x 'as 



can be simulated. Then 



[H 



w ■ 



H 



>}=H 



X 1 



H, 



w t M (e ^ + w t M ( - ei) 2/<Mi e<) 

where i = 1, . . . , d, Xi ^ j/j or Wi ^ Wi for all i. 

2. m — 2n even: d interactions of the above 

form where Xi — —yi,Wi ^ —Wi for all i and 2 d 

interactions of the form H^ i=1 CzC ^ , Ci G {0, ±1}. 



and H^r ■* can be obtained according to theorem[TJ Now 
we consider boxes with edge length bigger than 2 assum- 
ing that H [ x ' ) and * ] can be constructed for all 



L x 11 w 

v& G B z , and let v {z+1 ^ = (v\ 



Xz+l) 



)eB. 



z+l- 



First we show that there exist p^ z \q^ z ' G B z such that 
p( z )+q( z ) = u (*+i) )jJ W_grW e B z (see figure^). There- 



fore we define the set Z„ 



{i\\v[ z+1> \ = z + l}. If we 
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FIG. 2: Boxes B z and B z+1 . Any vector t/ z+1) £ B z+1 can 
be decomposed into a sum of vectors from B z , i.e., v^ z+1 ' = 
p (z) + where p (z \q (z) € B z , such that - p {z) e B z 
(see proofs of Thms. [2] and [4}. 



take = J2i e z v &\ e i and P {z) = v(z+1) ~ 1 (z) then b y 
definition pW,g-W G S z and p( z ) + gW = v( z+1 \ As for 
all i G Z v we have |(pW - q^). t \ = (H - 2) ^ z - 1 
and for all z ^ Z v we have |(p^ — g^),| = ^ z, it 
follows that p( z ' — £ B z . Using now the commutator 
relation 



i H x > H w \- l \ H w 



w 

and from (f2"6")) we know that we also get H { x 



(see Eq.([28])) we obtain that if, 



can be simulated, 

r(- (2 + 1) ) 

'x 

Now consider the case where m = In. If Xi — —y%, u>o ^ 
—Wi Vi, we can simulate and H^ e ^ for all i and 

A: according to theorem [1] (j2"c)) . Like in the case m = 

2n + 1 we can simulate H<£ v ] for all G Si, but 

according to theorem Q] (|2b[) simulating H^" ■* seems 
not to be possible. So we include these nearest-neighbor 
interactions in our initial set. The rest of the proof is 
then like in the case m — 2n + 1, and we see that the 
space 1Z can be simulated. □ 



and their commutator L" = [L, L'\ is given by 
L" = 



—W 
-X" 



X" 



Y" 
W" T 

X'(W + W T ) + X(W + W' T ), 



Y" = Y(W 



w" 



w 



IIT 



- W lT ) - Y'{W + W T ), 
X'Y - XY'. 



(30) 

(31) 
(32) 
(33) 



As in the fcrmionic case all commutators obey a symme- 
try which is in this case W" = W" T corresponding to re- 
flection symmetry (point symmetry) of the Hamiltonian 
and we will denote the vector space of all point symmetric 
Hamiltonians by V . This means that all simulated inter- 
actions are point symmetric up to linear combinations of 
the initial Hamiltonians. 

Every generator L of an arbitrary on-site Hamiltonian 
is of the form 



E, 



{xe,Ve,u)e) 



-we J Ue^ 



, x Ei Ve,w e G R. (34) 



Generators corresponding to a nearest-neighbor interac- 
tion along an axis are of the form 



-W Y 
-X W T 



(35) 



xM. 



(ei) 



where x, y, w, w € R and M'"' has been defined in JT 



We define L 



(v) 

w 



( 



Mi v) )®Mi v) and 





-m£° 



-(v) 



M. 




(v) 



(36) 



where the indices X, Y and W correspond to a non-zero 
X-, Y- and W-block respectively. 



A. Simulations in one-dimensional bosonic systems 

In this section we show that for one-dimensional 
bosonic systems an arbitrary nearest-neighbor interac- 
tion is complete for simulating the vector space V (i.e. 
W = W T ) when supplemented by all on-site Hamiltoni- 
ans. 



V. SIMULATIONS IN BOSONIC SYSTEMS 

In the following section we will study simulations in 
translationally invariant bosonic systems using quadratic 
on-site Hamiltonians and nearest-neighbor interactions. 
According to Eq. (fT3"|) the generators are of the form 



Theorem 3 Consider bosonic systems with quadratic 
Hamiltonians on a one- dimensional translationally in- 
variant lattice with periodic boundary conditions. The 
set of all possible one-mode transformations Ei XEtVEtWB \ 
in (|34p together with one arbitrary nearest-neighbor in- 
teraction given by L in Eq. (I35|) is universal for simulating 
the space V of all point symmetric interactions. 



L = o 1 



-W Y 
-X W 



Proof First we will show that an arbitrary interaction 
with X = Y = 0, W = W T can be brought from the 
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VF-block in the X and F-block: 



VI. SIMULATIONS IN SPIN SYSTEMS 



-W 

w 

-W 

w 





"5 



w 




-W 



(37) 
(38) 



Thus it is sufhcient to show that an arbitrary VF-block 
can be obtained. Let us start with a nearest-neighbor 



interaction of the form L 



(ei) 



As 



[Ly 1 ', £(1,0,0)] — L 



(ei) 



we also get L^ l} according to ([37]). Now [£^ lj ,£^ i; ] 



(ei) r( e i)l 



2£, 



(0,0,1) 



-(2ei) 



so that we also get L^ ei ^ 



and L 



(2ei) 

y 



As f£ 



(feei) r(ei)i 



y 



■((fe+l)ei) 



((fe-l)ei) 



we can simulate 



Finally it remains to show that we can get Ly 1 " 1 from 
an arbitrary nearest-neighbor interaction. If y =/= in 

then [[£ (ei) , £(0,0,1/2)], £(-i,o,o)] = so that 

we get Ly 1 ^ according to (|3"T|) . If y = 0,2: ^ 0, then 
[[£ (ei) , £(0,0,1/2)], £o,i,o] = xLfy' and we get Ly 1 ' as 
If s = y = 0,10 ^ 0, then [L^\ £ (0 ,_i,o)] = 



before 



^)4 ei) . 



□ 



B. Simulations in d-dimensional bosonic systems 

The following generalizes the previous result to cubic 
lattices in arbitrary spatial dimensions in cases where 
nearest neighbor interactions along all axes and diagonals 
are available. 

Theorem 4 Consider a system of bosonic modes on a d- 
dimensional translationally invariant lattice with periodic 
boundary conditions. The set of all on-site transforma- 
tions together with all nearest-neighbor interactions cor- 
responding to CiEi \ a S {0, ±1}, i = 1 . . . , d with L 
as in Eq. \35\) is complete for simulating the space V of 
all possible point symmetric interactions. 

Proof Like in the d-dimensional fermionic case let B z = 
{v|||«||oo < z}, z G N,v (z *> e B z . iFrom Thm. [5] we know 
that it is sufficient to show that L^/ can be simulated 
for arbitrary v. By induction over z we will show that 

i[y can be simulated for all S B z . For z — 1 
we know from Thm. [3] that all interactions described by 
G £i can be simulated as we have chosen our initial 

Hamiltonians appropriately. Now assume that ■* can 
be simulated for all t>w G B z , and let i>( z+1 ) G B z+ \. 
Then there exist p^ z \q^ G B z such that p^ + q( z > = 
v (z+i) )p (z) _ q (z) e Bz ( see n g ure [2]). As 



[L 



we can simulate L 



(p«) r (,W), = 



X 



J W 



J W ' 



ir 



□ 



In this section we will consider translationally invariant 
quantum lattice systems where a D-dimensional Hilbcrt 
space is assigned to each of the sites. We refer to these 
systems as spins although, of course, the described de- 
grees of freedom do not have to be spin-like. The main 
result of this section is that within the translationally in- 
variant setting universal sets of interactions cannot exist. 
These results are based on the following Lemma involving 
Casimir operators, i.e., operators which commute with 
every element of the Lie algebra [H, EH : 

Lemma 5 Consider a Lie- Algebra C and subalgebra 
£ = [C,£\. Let Sc be a set of generators for C and 
C a Casimir operator of C fulfilling 



ir[CG] = 0VG eS c \C. 
Then for every K G C we have that tr[CK] = 0. 



(39) 



Proof Every K G C can be written as 

K=Y,a L L+ Y] G G, a L ,/3 G GC. (40) 



Lea 



E 

Ges c \c 



Since we can write any L G £ as L = [£i,£a], Li G C 
we have that 

tv[CL] = tr[C[Li,£ 2 ]] = tr[[CL 1> L 2 \] = 

where we have used that C is a Casimir operator, i.e., 
VX G C : [C,L] = 0. Hence if we take the trace of 
Eq.gD]) with C we get 

tx[CK] = 0otr[CG] 
Ges c \c 

which vanishes according to the assumption in Eq. (f3"9")) . 

□ 

Let us now exploit Lemma [5] in the translationally in- 
variant setting in order to rule out the universality of 
interactions corresponding to certain sets of generators 
Sc- The following results are stated for onc-dimcnsional 
systems but they can be applied to d-dimensional lattices 
by grouping sites in d — 1 spatial dimensions. 

We use Casimir operators of the form 



k=l 



(41) 



where T\i l7 i 2 , ... ,i m ) = \h, *3, • • • , im, h) is the trans- 
lation operator which shifts the lattice by one site. To 
simplify notation we define an operator 



t(X) = Y T j XT^ 



(42) 



which symmetrizes any operator X with respect to the 
translation group. If X does not act on the entire lattice 
we will slightly abuse notation and write X instead of 
X ® 1® ••• ® 1. 
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Theorem 6 Consider a translationally invariant spin 
system on a ring of length to. If f is a non-trivial fac- 
tor of to then there is no universal set of Hamiltonians 
with interaction range smaller than f which generates all 
translationally invariant interactions. In particular if m 
is even, nearest-neighbor interactions cannot generate all 
next-to-nearest neighbor Hamiltonians. 



Now consider a two-body interaction between site one 



and site / + 1 of the form H = r{a - a 
Eq.dHJ we get 



tr[Off] = mii[af ] af +f) ] f[ tr[l] 

0=2 



From 



2mD f 



Proof Let us introduce a basis of the Lie algebra su o m 
of the form iaj 1 (8 ■ • ■®<Jj m , jk G %d 2 where o~k is traceless 
except for ao = 1 and tr[cr|] = 2 for all k > 0. For D = 2 
these are the Pauli matrices 



1 

1 



0-2 



-i 

1 



i 0- 3 



1 

-1 



(43) 



and for D > 2 we can simply choose all possible embed- 
dings thereof. Let m = /•/' and consider a Casimir oper- 
ator of the form in Eq. (|4"Tj) with 7^ = 6k, f, i.e., C =T* . 

We first show that for every Hamiltonian H G suo m 
with interaction range smaller than / we have ti[CH] = 
0. To see this note that the shift operator T* contracts 
the trace of a tensor product as 



tr[T^ 



ntr 

0=1 



f'-l 

a=0 



(44) 



In order to arrive at formula (|44[) expand the translation 
operator T* in the computational basis 

T f = ^2 \ i l--- i m)(i m -(f-l)irn-(f-2)---il---im-f\ 



in order to get 

tr [T f a jl <g> ••• ® a jm ] = 

( i m-(f-i)Wh\ii)(i m -(f-2)\o-j 2 \i2) ■ ■ ■ (i m -f\o- jm \i m ) 

i\ ,. . .,im 

Rearranging the order of the factors and using (/' — 1)/ = 
m — f leads to 

II ( i m~(f-0)\o-j ff \i0)(if3\o- ju+f>) \if+0) ■ 

0=1 ip,i f+p ,... 

(i(f>-2)f+0\o-j af ,_ 1)f+m \i(f-i)f+0) = 



ntr 

0=1 



Q = 



3(.0+<*f) 



If the interaction range of H is smaller than / then it 
can be decomposed into elements of the form r(ai 1 ® • • • ® 
a if). Since <7fc is traceless for all k > 0, all these terms 
amount to a vanishing trace in Eq.(|44p so that we have 
indeed tv[CH] = 0. This means we can apply Lemma [S] 
to the set Sc corresponding to all interactions with range 
smaller than /. 



such that by Lemma [5] we conclude that H cannot be 
simulated. □ 

The following shows that a universal set of nearest- 
neighbor interactions cannot exist irrespective of the fac- 
tors of to: 

Theorem 7 Consider a ring of length to. Then the 
set Sc corresponding to all on-site Hamiltonians and 
nearest-neighbor interactions is not universal for simu- 
lating all translationally invariant Hamiltonians. In par- 
ticular for D — 2 a product Hamiltonian 



H = t a 



(45) 



cannot be simulated if ct\, 02 and 03 all occur an odd num- 
ber of times. 

Proof We use the Casimir operator C = T — and 
the set of generators Sc = «{ T ( cr fe ® c su 2 ™. As 
tr [CG] = for all G G Sc we can again apply Lemma 
[5] Consider now the above product Hamiltonian H or if 
D > 2 its embedding respectively. Using Eq. (|4~4"|) with 
/ — 1, /' — m we obtain 



tv[CH] = mtr 



IK -IT 



k=l 



1 = 1 



(46) 



Since af — (— l) Si:2 ai and by assumption 01 appears an 
odd number of times we get tv[CH] = 2mtr [n)^ crjj 
which is non-zero iff 01 and 03 appear and odd number 
of times as well. □ 

Clearly, one can derive other no-go theorems in a simi- 
lar manner from Lemma O However, we end this section 
by providing some examples of interactions which can be 
simulated. For this we define gki — r(ak <8> 07). 

Theorem 8 Consider a translationally invariant system 
of m qubits (D — 2) on a ring. By using on-site Hamil- 
tonians and nearest-neighbor interactions the following 
interactions can be simulated: 



4 3) = 



1 0~i 
1 (To 



(47) 
(48) 



where i,j G {1,2,3} and rj denotes the number of aj 
matrices. Moreover, for m = 5 one can simulate next- 
to-nearest neighbor interactions of the form N{ = r(crj (g) 
1® Oi ® 1(g) 1). 
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Proof We will restrict our proof to the pairs i = 1, j = 2 
as the other interactions can be obtained in an analo- 
gous way. We start proving (|47jl . The Hamiltonian |47|) 
can be simulated due to [321,513]/ (2«) = r(cr 1 (g) (g ci ) . 
For proving (|^5|) . we start with [313, gu]/(2i) — . As 
[J^\g 31 ]/(2i) = - J[f = 9ll , we have 

shown (g5]). 

Now we will prove that the next-to-nearest neighbor in- 
teraction can be achieved for m = 5, i = 1 (i = 2, 3 follow 
similarly) . Using ([48]) , we see that t (01 <g eri <g ci <g en <g 1) 

can be simulated, as [323, J21 ]/(2i) + ^(21) — T (°i ® °"i ® 
cti <g ci <g 1) , and similarly we get r (0-3 03 <g 03 (8 0-3 <g 1) . 
As [[311,323], 532] /4 = 2Ni - r{a l <g a x ® o x ® a x <g 1) - 
r((73 (g 03 CT3 <g 03 ® 1) we can extract N±. □ 



VII. CONCLUSIONS 

We have presented a characterization of universal sets 
of translationally invariant Hamiltonians for the simula- 
tion of interactions in quadratic fermionic and bosonic 
systems given the ability of engineering local and nearest 
neighbor interactions. Thereby the Lie algebraic tech- 
niques of quantum simulation restrict the space of reach- 
able interactions to Hamiltonians with real hopping am- 
plitudes in the case of fermions and to point symmetric 



interactions in the case of bosons. 

For spins the situation appears to be more difficult and 
a complete characterization of interactions that can be 
simulated remains to be found. As a first step, we have 
identified Hamiltonians that cannot be simulated using 
short range interactions only. Furthermore, we have in- 
troduced a technique based on the Casimir operator of 
the corresponding Lie algebra which allows one to find 
Hamiltonians that cannot be simulated with a given set 
of interactions. 

In this work we have considered the question of what 
can be simulated leaving aside the question of the effi- 
ciency. In this context it is important to remark the fact 
that the number of applications of the original Hamilto- 
nians in order to obtain a result bounded by some given 
error scales polynomially in the Trotter expansion. The 
scaling with the total number of particles depends on the 
number of commutators that are required to obtain the 
Hamiltonian. 

Finally, whereas we have shown that it is not possible 
to perform certain simulations for spin systems it is still 
possible to perform those simulations by encoding the 
qubits in a different way. 

We thank the Elite Network of Bavaria QCCC, DFG- 
Forschungsgruppe 635, SFB 631, SCALA and CON- 
QUEST. 
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